Approximating chiral quark models with 
linear <r-models* 



Wojciech Broniowski a and Bojan Golli b 

a The H. Niewodniczanski Institute of Nuclear Physics, ul. Radzikowskiego 152, 

PL-31342 Krakow, Poland 

b Faculty of Education, University of Ljubljana and J. Stefan Institute, Jamova 39, 
P.O. Box 3000, 1001 Ljubljana, Slovenia 



Abstract 

We study the approximation of chiral quark models with simpler models, obtained 
via gradient expansion. The resulting Lagrangian of the type of the linear u-model 
contains, at the lowest level of the gradient-expanded meson action, an additional 
term of the form ^A(ad fl a + ird^ir) 2 . We investigate the dynamical consequences 
of this term and its relevance to the phenomenology of the soliton models of the 
nucleon. It is found that the inclusion of the new term allows for a more efficient 
approximation of the underlying quark theory, especially in those cases where dy- 
namics allows for a large deviation of the chiral fields from the chiral circle, such as 
in quark models with non-local regulators. This is of practical importance, since the 
cr-models with valence quarks only are technically much easier to treat and simpler 
to solve than the quark models with the full-fledged Dirac sea. 
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1 Introduction 



Various cr-models, i.e. effective models incorporating the spontaneously broken 
chiral symmetry, play an important role in modeling the low-energy hadronic 
physics. The original Gell-Mann-Levy linear cr-model [1], after replacing nu- 
cleons by quarks [2,3,4,5,6,7], has been used to describe the binding mecha- 
nism of quarks inside baryons and to produce quite successful phenomenolo- 
gies of baryons. Actually, these models, with Lagrangians containing up to 
two derivatives of fields, may be viewed as formal approximations to models 
with purely quark degrees of freedom, such as the Nambu-Jona-Lasinio (NJL) 
model [8] with quarks (for review and further references to these models in 
the context of mesons see, e.g., [9,10,11], and in the contexts of chiral soli- 
tons e.g., [12,13,14,15]). Indeed, gradient expansion techniques [16,17,18,19] 
produce mesonic terms of the form used in various cr-models. Admittedly, 
questions of convergence of the gradient-expansion series are difficult, and it 
is not clear how many terms should be present in the mesonic theory such that 
it approximates satisfactorily the underlying quark theory. On the other hand, 
the cr-like models are much simpler to solve, since they do not carry the com- 
plications of the Dirac sea, which for models of baryons is technically difficult 
to treat, and a massive computational effort is involved. Therefore, it is desir- 
able to understand how to construct much simpler and more intuitive effective 
a-models which carry valence quarks only, and then to investigate what role 
is played by various mesonic terms arising from the gradient expansion. 

In this paper we analyze a particular extension of the linear cr-model coupled 
to valence quarks, which contains an additional term with two gradients of the 
chiral fields. This term, arising naturally in the process of gradient-expanding 
NJL-like models, has not been, to our knowledge, considered before. It has 
dynamical consequences for the model and, in addition, allows for a more 
efficient approximation of the underlying quark model. We discuss the nature 
of the term and its consequences for the Euler-Lagrange equations in Sec. 2. In 
Sec. 3 we show how the term arises when the NJL model is gradient-expanded 
and explain how the parameters of the resulting cr-model can be derived. We 
estimate the strength of the new term in NJL model with the sharp four- 
momentum cut-off, and find that it is large. We then proceed to analyzing its 
dynamical effects for chiral solitons. In Sec. 4 we consider the straightforward 
extension of the cr-model, and find that the new term lowers the energy of 
the soliton, thus contributing to binding and its energetic stability. In Sec. 
5 we apply our analysis to the chiral quark model with non-local regulator 
[14,15,20,21,22,23,24,25,26,27,28], where we show that the inclusion of the 
new term allows for a more efficient approximation of the underlying quark 
theory. 
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2 The A-term 



By means of general symmetry arguments, in the construction of the effective 
Lagrangian one has at hand the invariants of the given symmetry group, in 
the case considered here the chiral SU{2)l x SU(2)r. The chiral fields, i.e. the 
scalar-isoscalar cx-field and the three pion fields, ir a , belong to the (1/2, 1/2) 
representation of this group. The typical chiral invariants one can form out of 
the chiral fields, arranged in the increasing number of spatial derivatives, are 

a + tc , 

ad^a + 7T a d^n a , 

(d»a)(d v a) + (d^7i a )(d u 7T a ), 

(8 fia a)(d ,/ (T) + (d> ia n a )(d u TT a ), (1) 
etc. 

Similarly, the coupling to Dirac fields, if), allows to construct invariants of the 
form ip(<j + i r )^T a "K a )%l) 1 ^7 M (o"9 M o" + TT a d' J "ir a )ip, etc. The above invariants are 
next combined in such a way as to form Lorentz invariants. In the mesonic 
sector, up to terms involving two gradients, we have the general form for the 
mesonic Lagrangian, 

L mes — —V + \z[{d^f + (c>X) 2 ] + \A[ad»cj + n a dV] 2 , (2) 

with V, Z, and A denoting, at the moment, arbitrary functions of the invariant 
cr 2 + 7r 2 . The last term in Eq. (2), called from now on the A-term, is new 
compared to previous works. In the original Gell-Mann-Levy model this term 
is excluded by the requirement of renormalizability. Since we are going to 
use (2) as an effective model, approximating the underlying quark theory, the 
model need not and should not be renormalizable, hence the A-term may be 
present. Moreover, as we will show in Sec. 3, the appearance of the A-term 
naturally follows from the gradient expansion of the quark models with the 
Dirac sea, such as the NJL model and its modifications. 

The linearity of the model is an essential feature in our study. The non-linear 
constraint imposed on the chiral fields, a 2 + tt 2 = const., immediately kills the 
A-term, since ad^a + 7i a d' l 7i a = ~ ^((j 2 + n 2 ) = 0. As a matter of fact, the 
underlying quark models lead to linear chiral models, and the non-linearity in 
most approaches [12,13,14,15] is superimposed externally. Q 

1 In fact, in these works the non- linear constraint is needed to provide stability of 
solitons [12,13]. On the other hand, in non-local models of Refs. [24,25] stability of 
solitons is achieved without the "external" non- linear constraint. 
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In the following we shall use the coupling to quarks of the form as in the Gell- 
Mann-Levy model [1], i.e. —gip(a + i r )^ ) T a 'K a )ip ) with ip containing the valence 
quark orbit only. The Dirac sea has been (approximately) integrated out and 
entered into the mesonic part of the Lagrangian. In addition, for simplicity of 
the calculations, we assume that the functionals Z and A are constant, i.e. do 
not depend on the value of the chiral fields: 



Z(a 2 + TT 2 ) = Z , A{a 2 



7T 



Aq. 



(3) 



This assumptions is justified as long as the field do not depart too far away 
from the chiral circle. 

The mean-field Euler-Lagrange equations resulting from Eq. (2) have the form 



dV 

~do~ 

dV 



Ja = 0, 



Z Oa + A a [{d^af + (<V a ) 2 + ffDff + ^ u ^ a 
Z nn a + A n a [(c» 2 + (d,n b ) 2 + aUa + n b Dn b ] +^-+j% = Q, 
where j a and j° are the valence quark sources, 

idval idval 



(4) 



(5) 



and qi are the valence orbitals.Q As can be seen from Eqs. (4), the A-term 
mixes the propagation of the chiral fields. Equations (4) can be diagonalized 
and brought to the simple matrix form 



Uo | 1 

n7l a J = ~ Z (Z + + Ao) * (6) 

Z + A 7T 2 -A an a \ ( A Q a((d^a) 2 + (<V) 2 ) + dV/da + j a 
-A a7T a Z + A a 2 ) { A -n a {{d^af + (<V) 2 ) + dV/dn a + # 

If A Q = 0, the usual Gell-Mann-Levy a-model equations of motion are recov- 
ered. 

Equations (4) are simple to solve in the vacuum sector. In this case the A and 
Z terms do not contribute and one has the familiar result: 

(a) = F, (7T a > = 0, (7) 



The effects of the Dirac sea are included in the mesonic degrees of freedom. 
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where F = 93 MeV denotes the pion decay constant. As usual, we introduce 
the shifted u-field, a' = a — F, and find for the small-amplitude fluctuations 
of the chiral fields around their vacuum values (7) the following equations of 
motion 



(Z + A F 2 )Da' + 



d 2 V 



da 2 



a' = 0, 



=7T=0 



Z an a + 



d 2 V 



dn a dn h 



TV 



0. 



(8) 



= 7T = 



We can thus canonically normalize the pion field by choosing Zq — 1. Then, 
for non-zero Aq, the normalization of the a field is not canonical, and carries 
the factor 1 + A F 2 , which means that the mass squared of the physical a' 
field is given by the expression 



;i + A F 2 ) 



da 2 



(9) 



=7T=0 



One can easily check that the A-term does not contribute explicitly to the 
vector, V£, and axial-vector, A a , Noether currents. The meson contributions 
remain to be the usual cr-model expressions, 



V « = Ze abc (d^ b )7r c , 

A; = Z(ad,7T a -n a d^a). (10) 

The above expressions hold for the case of general Z and A. Therefore the 
expressions for the magnetic moments or qa used in soliton calculations [29] 
are not modified. However, changes in these quantities are induced by the 
A-term through the change of the dynamics in Eqs. (6). 



3 Motivation from quark models 



We are now going to argue that the presence of the A-term in effective chiral 
Lagrangians follows naturally if one considers approximations to chiral quark 
models, such as the NJL model. For simplicity of notation we work in the strict 
chiral limit of the vanishing current quark mass. The action of the bosonized 
NJL model at the one-quark-loop level can be written as (see, e.g., [12,13,15]) 

I(o-,n a ) = -Tr\og(d T + h)-^Jd*x(v 2 + * 2 ), (11) 
h = -ia- V + pg(a + i-f 5 r a n a ). (12) 
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The trace is over color, flavor, Dirac, and functional space, r denotes the 
Euclidean time. We have introduced an additional parameter, g, the quark- 
meson coupling constant, which allows us to fix the vacuum value of the a- 
field by imposing the condition (7). The constant G is the four-quark coupling 
constant of the original N JL model [8] . The constituent quark mass is related 
to the pion decay constant, 

M = gF. (13) 



For many practical applications (study of large-size solitons, checking numer- 
ical computations) it is useful to have an approximation to this model, which 
is much simpler to deal with. In practice, this amounts to the gradient expan- 
sion [16,17,18,19] of action (11), with the hope that for sufficiently large-size 
solitons the expansion should work fine. Also, in many models the fields lie 
close to the chiral circle, therefore another good expansion parameter is the 
combination 

5 = a 2 + 7t 2 - F 2 , (14) 



which measures the departure of the fields from the chiral circle. We can thus 
expand 

oo oo oo 

V(6) = £ VJ n , Z{5) = £ ZJ n , A{5) = £ AJ n . (15) 

n=0 n=0 n=0 



The coefficients V n , Z n and A n can be found by comparing the Green's func- 
tions obtained from (11) to those obtained at the effective level from (2). 

Our technique may be viewed as a double expansion: in the number of gra- 
dients, up to two, and in the parameter 5. In the functions Z and A we thus 
keep the terms up to zeroth order in 5, while in the potential V we shall keep 
term up to second order in 5. 

To begin, we recall the gradient expansion of V, which is easy, since for that 
purpose one can treat the fields as space-time independent. We also have the 
grand-reversal symmetry [30], which leads to Tr log (d T + h) — Tr log (— d T + h) = 
|Tr log (-d 2 + h 2 ). Therefore 



V(a 2 + n 2 ) = -2N c N f J ^ log (k 2 + g 2 (a 2 + n 2 )) + ^(a 2 + vr 2 ) 

+ const., (16) 

where k is the Euclidean momentum flowing along the quark loop, iV c = 3 is 
the number of colors, and Nf = 2 is the number of flavors. The constant can 
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be chosen such that V = in the vacuum (7). In addition, the stationary-point 
condition, dV/da\ Yac = 0, has the explicit form 

1 » (17) 



2G 2 J (2vr) 4 A; 2 + M 2 



Plugging (17) into (16) results in the following formula for V: 



V(5) = -2N c N f J 



d 4 k 



g 2 5 



log 1 + 



9 2 S 



k 2 + M 2 k 2 + M 2 



(18) 



The coefficients of the expansion (18) are: 



We should bear in mind that the expansion for the log is rather slowly con- 
vergent, such that many terms in Eq. (15) may be needed to reproduce well 
the full result (16). Retaining the n = 2 term only corresponds to the use of 
the usual Mexican hat potential. 

Next, we calculate the coefficients Z and A . For that purpose one can con- 
sider mesonic two-point Green functions. In the vacuum the inverse a and 
pion propagators obtained from Eq. (2) are (the momentum q is Euclidean): 



5 2 L 



5n a (q)5'K b (-q) 

5 2 L d 2 V(a 2 + 7T 2 

5a(q)5a(-q) 



Z q 2 5 ab , (19) 
+ (Z + F 2 A ) q 2 . 



da 2 



Meanwhile, in the NJL model one obtains, with help of Eqs. (11,12), the 
following formulas for the inverse meson propagators: 



K^ = 2g 2 N c N f f(M,q 2 )q 2 = 2g 2 N c N f f(M, 0)q 2 + 0(q 4 ), (20) 
K; 1 = 2g 2 N c N f f(M,q 2 )(q 2 + 4M 2 ) = 8g 2 N c N f M 2 f(M, 0) 

+ 2g 2 N c N f (>(M,0)+4M 2 df(M,q 2 )/dq 2 \ q= ^q 2 + 0(q 4 ), 

where 

W, q 2 ) = I — 9 — 9 _ v (21) 



(2vr) 4 ((jfe + q/2) 2 + M 2 ) {[k - q/2) 2 + M 2 ) 
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The index A in the integral denotes an appropriate regulator [8], necessary 
for the model to become finite. By comparing (19) and (20) we immediately 
arrive at the relations 



Z = 2g 2 N c N f f(M,0), 

A = 8g 2 N c N f (M/F) 2 df(M,q 2 )/dq 2 



(22) 



g=0 



Through the use of the relation F = M^2N c Nff(M, 0) [12] and Eq. (13) we 
obtain the formulas 



Z — l or g~ 



2N c N f f(M,0), 



(23) 



, , 2 df(M,q 2 )/dq 2 \ q=0 
An = 4q ; — - — . 



In order to estimate the size of the A-term, in Fig. 1 we plot the quantity 
1 + F 2 A , which plays the role of the wave function renormalization for the 
a-field. In the study of this section we have used the sharp four-momentum 
cut-off on the quark loop in Eq. (21). The calculation is performed in the 
chiral limit, and with the cut-off A chosen in such a way as to reproduce the 
physical value of F = 93 MeV. We note that the effect of the A-term is very 
strong, with the combination 1 + F 2 An ranging from 0.6 at M = 300 MeV to 
as little as 0.1 at M = 600 MeV, to be compared to unity of the case without 
the A-term. 

It is easy to derive full expressions for Z and A. With the notation of Eq. (14) 
we find 



Z(S) = 2g 2 N c N f J — -i— — ( 24 ) 
A(5) = -2g*N c N f I A^/ + H^ + M 2 ) 

^ ' y i J (27r) 4 (fc2 + ( 5 2 tf + M 2))4 V ) 

These formulas can be straightforwardly expanded in powers of 5 to produce 
all coefficients in (15). 

Note that while the sign of Z is positive-definite, the sign A is negative- 
definite. Therefore the presence of the A-term in the effective Lagrangian 
lowers the wave function renormalization of the <T-field, thereby increasing its 
mass, as follows from Eq. (9). 
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4 Chiral solitons 



In the following parts of this paper we are going to explore the dynamical 
consequences of the A-term for chiral solitons. First, we consider the simplest 
case: the generalization of the Gell-Mann-Levy Lagrangian of the form 

L = \[{d,af + (c> a ) 2 ] + ~A [ad»a + 7r a dV} 2 - ^(a 2 + vr 2 - v 2 ) 2 
+ ca- g^p{i$ + a + i7 5 tV)V>, (26) 

with 



— 2 2 

mi — mi 



mi 



2F 2 



v 



mlF, 



mi 



m 2 



[l + F 2 A )ml. 



(27) 



The last equality ensures that the sigma mass is m a , according to Eq. (9). 
The constants g, m a , and A are treated as model parameters. The resulting 
mean-field equations of the form (6), and the equations for the upper and lower 
components of the valence quark orbital [2,3], are solved with the hedgehog 
ansatz for the chiral fields: a(x) = cr(r), ir a (x) = r a jr 7r(r), where r denotes 
the radial coordinate. 

A typical solution, for g = 4.5, m a = 900 MeV, and F 2 Aq = —0.6 (dashed 
lines) and Aq = (solid lines), is presented in Fig. 2, where we show the 

1 

0.8 

CM 

LL 

<°0.6 
+ 

0.2 

300 350 400 450 500 550 600 

M [MeV] 

Fig. 1. The quantity 1 + AqF 2 evaluated in the Nambu-Jona-Lasinio model with the 
sharp four-momentum regulator, plotted as a function of the quark mass, M. Large 
deviation from unity is found. The calculation is performed in the chiral limit, and 
with the cut-off chosen in such a way as to reproduce the physical value of the pion 
decay constant. 
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chiral fields, and the radial valence quark density (radial baryon density), 
47rr 2 ps( r ) — 4:7rr 2 q^q, plotted as functions of the radial coordinate r. The 
chosen value of A is suggested by the analysis of Sec. 3, cf. Fig. 1. The 
values of g and m a are the typical values used in other studies. We notice that 
the presence of the A-term weakens the pion field by about 20%, modifies 
the shape of the a field, as well as slightly increases the size of the soliton, 
which is visible from the baryon-density curves. The results do not depend 
qualitatively on the values of the parameters g and m a . As expected, the 
effect is more pronounced if the chiral fields lie away from the chiral circle. 

In Fig. 3 we study the behavior of the soliton energy as a function of the 
coupling constant g. Again, we use F 2 A = —0.6 (dashed lines) and A = 
(solid lines). For m a we take two "extreme values", 400 MeV and 1200 MeV. 
We note that the inclusion of the A-term lowers the energy of the soliton 
for all values of g. The lowering is a significant effect, by the amount of a 
few tens of MeV. The solitons become energetically stable when their energy 
is lower than the energy of three free quarks, 3M, denoted in Fig. 3 by the 
straight solid line. With the A-term present, the stability is achieved at lower 
values of g. This is a desired effect since many phenomenological approaches 
[31,32,33,34,35,36] have problems in getting the mass in the right ball park if 
g is too high. 

— i — i — i — i — i — i — i — i — i — | — i — i — i — i — i — i — i — i — i — | — i — i — i — i — i — i — i — i — i — 




Fig. 2. The chiral fields and the radial quark density in the hedgehog chiral soliton 
of model (26), plotted as functions of the radial coordinate r. Two solutions are 
compared: with Aq = (solid lines) and AqF 2 = —0.6 (dashed lines). The remaining 
model parameters are g = 4.5 and m a = 900 MeV. 
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On the other hand, as expected, the A-term has negligible effect on the nucleon 
observables (only of the order of few percent), since it does not explicitly enter 
in the expressions for the Noether currents (10). 



5 Model with the non-local regulator 



Recently, considerable progress has been made in the treatment of chiral quark 
models with non-local interactions. It these models, rather then cutting the 
momentum in the quark loop, the finiteness of the theory is achieved by the 
non- locality of the four-quark interaction. In fact, the model derivations of 
effective chiral quark models, such as the instanton-liquid model [14,37], or the 
resummations of rainbow diagrams in the Schwinger-Dyson approach [38,39], 
lead in a natural way to non-local models [15]. Their applications range from 
the meson sector [22,23,40,41], through quark matter at non-zero temperature 
[42] and baryon number [26,43,44,45], chiral solitons [24,25], to applications 
in exclusive processes in QCD [27,28]. The advantages of non-local models 
include such features as the natural preservation of anomalies, finiteness to all 
orders in the l/iV c -expansion, such that the meson-loop calculations can be 
done with no extra parameters [41,42,46], or the stability of solitons without 
the need for the non-linear constraint. The solitons of non-local models, their 




1400 - 



800 - 



_l I I I I I I I I I I I I I I I I I I I I I I l_ 



3.5 



4.0 



4.5 



5.0 



Fig. 3. The dependence of the energy of the soliton on the coupling constant g for 
Aq = (solid lines) and AqF 2 = —0.6 (dashed lines). Two values of m a are used. 
The straight solid line is the energy of three free quarks. The presence of the A-texm 
noticeably lowers the energy of the soliton. 
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construction and the resulting phenomenology, have been discussed in detail 
in Refs. [24,25]. Here we are going to check how well we can approximate the 
full, complicated model, with the cr-type model supplied with the A-term. 



The effective action for the non-local model has the form (11), with the crucial 
difference that now the Dirac Hamiltonian carries the regularization operator, 
r [15]: 



h = -ia ■ V + r(d 2 )(3g(a + i lh T a n a )r(d 2 ) + (3m. 



(28) 



The current quark mass is denoted by m. It is convenient to introduce the 
notation R = r 2 . One can straightforwardly derive the expression for the 
effective potential, 



V(a,ir a ) = -2N c N f J 



x 



log 1 + 



d±k 

g 2 R 2 5 + 2mgR(a - F)\ g 2 (R 2 + Rm /M) 8 



(29) 



k 2 + (RM + my 
1 



k 2 + (RM + mf _ 



V 2 (a 2 + n 2 -F 2 ) 2 + ^mi (a — F) + ir + 0(m\mb\ 5"), 



which gives 



V 2 = iV c iV^ 4 / 

jy J (2ir) 4 (k 2 + R 2 M 2 ) 



R 4 



(30) 



The formulas for Z and A obtained from the non-local model read 



2N. 



Nf9 2 J 



d^k R 2 - k 2 RR' + (k 2 R r ) 
(2tt) 4 (k 2 + R 2 M 2 f 



(31) 



Ao = -8N c N f g A J 



dik 

x 



(2nY 

R 2 ((k 2 + 2R 2 M 2 )R 2 - 2k 2 (2k 2 + R 2 M 2 )RR' + 8k 6 R' 2 ) 
4 (k 2 + R 2 M 2 ) 4 ' 



(32) 



where R' = dR/d(k 2 ). The model parameters are fitted in such a way that 
Zq = 1, which ensures that F has its physical value [22,23,25]. 

The presence of the regulator r dependent on the momentum modifies the 
Noether currents of the theory [24,25,47]. As a result, discussed in detail in 
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Refs. [24,25], the contributions of the valence quarks to the Euler-Lagrange 
equations carry an extra residue factor, z va x. This factor is defined as 



l-i 



dui 



(33) 



<^=*<:val, 



with e va i(o;) denoting the energy- dependent valence-orbital eigenvalue of the 
energy-dependent Dirac Hamiltonian (28). It satisfies the equation 



h (-<4i) l val ) = e vai|val). 



(34) 



Thus, 



2c (x) = gN c z v3l (val \r\x)/3(x\r\ val) , 

ft (x) = gN c z val (val| r | x)j3i^ b r a (x\ r \ val) . 



(35) 
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Fig. 4. The dependence of the energy of the soliton of the model with non-local 
regulators on the coupling constant g. The solid line on the right corresponds to 
the solution of the quark model with non-local regulators, the dotted line shows the 
solution of the approximating <r-like model with Aq = 0, and the dashed line shows 
the solution of the approximating model with the ^4-term present. The solid line on 
the left corresponds to the energy of three free quarks in the non-local model. 
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In the numerical study presented in this paper we use the Gaussian regulator, 
R(k)=r(ky = exp(-^- 2 ). (36) 

In Fig. 4 we compare the hedgehog soliton solutions of the non-local quark 
model, defined by Eqs. (11,28,36), with those of the approximating a-like 
model, defined by Eqs. (2,30,32,36), and prescriptions (33,34,36) for the va- 
lence orbit. We plot the energy of the solution as a function of the quark-meson 
coupling constant, g. The solid line on the right corresponds to the non-local 
quark model, the dotted line to the approximating model with A = 0, and 
the dashed line to the approximating model with A given by Eq. (32). Sim- 
ilarly to the study of Sec. 4, the inclusion of the A-term results in lowering 
of the energy. We note that the model with the A-term approximates the full 
non-local quark model better than the model with A = 0. In particular, at 
lower values of g the agreement is very good, since in that region the solitons 
are large and the gradient expansion works well even at the lowest order. The 
solid line on the left corresponds to the energy of three free quarks in the 
non-local quark model [24,25]. 

Comparing the nucleon observables calculated in the approximating model to 
those calculated in the full model [25] we find that the valence contribution 
to g\ and the magnetic moment is reproduced within a few percents for g 
below 3.5 (corresponding to M below 325 MeV), and within 10 % for g > 3.5. 
The meson contribution to the magnetic moment is again reproduced well for 
g ~ 3.3, and is typically 20 % larger than the sea contribution in the quark 
model for g ~ 3.7, and more than 30 % for g > 5. On the other hand, the meson 
contribution to g& is a factor of 2 too large compared to the sea contribution 
already for g slightly above the threshold value g = 3; this factor rises to 3 at 
g ~ 5. The results do not change much if we switch-off the A-term; the only 
noticeable difference is in the meson contribution to the magnetic moment at 
larger values of g where the inclusion of the A-term improves the results by 
5 % - 10 %. This is a dynamical effect reflecting the tendency of the A-term 
to reduce the strength of the pion field, as can be seen already in Fig. 2. 



6 Conclusion 

We have analyzed the models approximating the chiral quark models with 
four-quark interaction. The models are constructed with help of the gradient- 
expansion technique, good for large solitons. We have found that the inclusion 
of the A-term, appearing naturally in the gradient expansion of the underlying 
theory, improves the approximating model. The term results in a significant 
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(a few tens of MeV) lowering of the energy of the soliton. This effect is par- 
ticularly important in the non-local model, where the A-term considerably 
improves the approximating model such that for lower values of the quark- 
meson coupling the energy follows almost exactly the energy obtained in the 
full quark model. We wish to stress that the variants of the a-model, which 
result from the approximation of the full quark model via the gradient ex- 
pansion are much simpler to solve than the original theory. They avoid the 
complications of the numerical treatment of the Dirac sea, and thus are useful 
and more practical in the studies of soliton models of baryons. 

The authors wish to thank Enrique Ruiz Arriola and George Ripka for several 
useful discussions. One of us (WB) wishes to thank the organizers of the 
Bled Workshops in Physics, where this research was initiated, for their kind 
hospitality and stimulating atmosphere. 
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